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ABSTRACT 

We report our search for optically thick radio cores in sixteen narrow-waist 
bipolar nebulae. Optically thick cores are a characteristic signature of coUimated 
ionized winds. Eleven northern nebulae were observed with the Very Large Array 
(VLA) at 1.3 cm and 0.7 cm, and five southern nebulae were observed with 
the Australia Telescope Compact Array (ATCA) at 6 cm and 3.6 cm. Two 
northern objects, 19W32 and M 1-91, and three southern objects. He 2-25, He 2- 
84 and Mz 3, were found to exhibit a compact radio core with a rising spectrum 
consistent with an ionized jet. Such jets have been seen in M 2-9 and may be 
responsible for shaping bipolar structure in planetary nebulae. 



Subject headings: planetary nebulae: general — stars: AGB and post-AGB — 
stars: winds, outflows — stars: mass loss — stars: evolution — radio continuum: 
stars 
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Introduction 



Planetary nebulae (PNs) show a remarkably broad range of often intricate morphologies. 
Butterfly nebulae, cla ssified as bipolar nebulae with a narrow pinched waist at the center 
( iBalick fc Frankll2002l ). are probably the most spectacular, and yet at the sarae tini e most 
puzzling. Under the classical interacting stellar winds (ISW) model (iKwok et al.lll978l ). a FN 
is formed when the fast ionized wind from its central star overtakes the material ejected earlier 
by the slow neutral wind and compresses the material into a dense shell. This produces a 
round PN if both winds are isotropic. A dense equatorial torus, which blocks the fast wind in 
that direction, is added to the interacting winds system to exp l ain how asym metric structure 
in PNs (e.g., bipolar) can be developed (IKahn fc Westlll985l : iBalickl 119871 ). This, however, 
can only produce bipolar nebulae with wide waists, not those with narrow "pinched" waists. 
It ha s been suggested that the cr eation of the latter ri ebulae requires a collimated fast wind 



[e.g., iSoker fc RappaportI 120001 ). ISoker fc Rappaportr s model invokes a binary star system 
to produce a collimated fast wind. The system consists of a mass-losing AGB star and a 
main-sequence or white dwarf companion. The AGB wind captured by the companion forms 
an accretion disk around the companion that loads and drives a collimated fast ionized wind. 

The existence of collimated outflows (COFs) has been suggested by numerous obser- 
vations, usually from the ionized gas images. For example, the highly collimated lobes of 
some PNs and PPNs (proto-plane tary nebulae) imply that their morphology is shaped by a 
collimated outflow (see review by lKwokll2003al . and references ther ein). The first di s cover y 
of strongly collimated high- velocity bipolar outflows in a PN was by lGieseking et al.l (119851 ) . 
in NGC 2392, the Eskimo nebula. They analyzed from spectroscopic observations of this 
nebula, particularly in the [N II] A6583 emission line, they found a set of bipolar jets with a 
expansion velocity of ~ 200 km s~^ and a small angle of collimation less than 10°. Although 
unexpected at first, the pervasiveness of COFs in PNs is now widely recognized. COFs have 
diverse characteristics. They are prominent in emission lines such as [N II] and [S II] and 
are observed as bipolar, multipolar, point-symmetric and jet-like outflows. The expansion 
veloc i ties o f these COFs range from several tens to a thousand km s~^ (see reviews of [Lopez 
20021 . l2003h . 



In addition to ionized gas, COFs have also been seen in neutral gas. Sever al authors have 
made detailed studies on CO emission lines of individual PNs and PPNs (e.g.. lForveille et al. 



19981 : iHuggins et al.ll2000l ). The ve locities for these CO out flows are found to be much smaller 
than those of the ionized COFs. iBujarrabal et al.l (120011 ) studied the CO emission from a 
large sample of PPNs as well as two yellow hypergiants and one young PN. They reported 
collimated fast outflows in 28 of 32 objects, indicating that such outflows are a common 
phenomenon in the early PN stage. 
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In the model proposed by lSoker fc RappaportI (I2OOOI ). the coUimated winds responsible 
for shaping butterfly nebulae are formed and coUimated close to the central star system. 
They are difficult to observe in optical wavelengths directly, because dust around the equator 
limits the depth into the center that optical observations can probe. The dust is optically 
thin to radio emission, which can therefore probe deep int o the center t o look for c entral 
jets, if the jets are at least partially ionized. Using the VLA. lLim fc Kwokl (l200Cll . l2003l ) have 
indeed found a coUimated wind in the form of a pair of ionized jet in the core region of the 
narrow-waist bipolar nebula M 2-9. 

M 2-9 has a bright central core and two cigar-shaped lobes extending from opposite 
sides of the core ~ 20" to the north and south. The lobes of the nebula contain three 
pairs of symmetric knots or bright features: two of them are with mirror symmetry and 
one pair with point/mirror symmetry. Because of its spectacular morphology, M 2-9 is 
one of the most well-studied planetary nebulae. Since its discovery in 1947, the core and 
the outline of the lobes have remained stationar y, but the knots and bright features have 



move d laterally from the west to the east (e.g., lAUen fc Swingsl Il972l : iKohoutek fc Surdej 



19801). Various models h ave been propos e d to e xplain the lateral motions of bright features 
of M 2-9. For example, lAUen fc Swingsl (119721 ) proposed that the changes were the result 
of a rotating beam of i onizing radiation fro r n the star shining through a matched pair of 
holes in orbiting clouds. iKohoutek fc Surdejl (Il980l ) noted that it is unlikely such model can 
result in positional changes only in the east-w est direction, an d instead suggested that the 
nebula was rotating as a whole. However, sinc e iGoodrichI (119911) conclude d that the distance 
of M 2-9 cannot be only 50 pc as proposed by IKohoutek fc Surdejl (Il980l ). the unreasonable 
large velocities of the knots have made their model problematic. 



Recently in explaining the structural changes of M 2-9. lDoyle et al.l (120001 ) and lLivio fc Soker 



(120011 ) both invoke a symbiotic star system to produce a mirror-symmetric jet. Direct evi- 
dence of such a jet i s provided by observations of M 2-9 at 6 cm, 3.6 cm, 2 cm, and 1.3 cm by 
(iLim &: Kwokl l2003l ). Their observations were largely insensitive to the large-scale free- free 
em ission from the neb ula and focused instead on the compact radio core at the center found 
by iKwok et al.l (119851 ). Their results indicate that this optically thick core has a spectral 
index of 0.67 ± 0.01, which is close to the spectral index of 0. 6 expected frorn an isothermal 
outflow expanding at a constant velocity and opening angle (iReynoldd Il986l ) . More impor- 
tantly, at the shortest wavelength of 1.3 cm wher e the jet is best res olved, its bipolar lobes 
show a mirror-symmetric structure (Figure 1 of iLim fc Kwokl l2003l ). Both lobes point in 
the directions of mirror-symmetric brightenings in the n ebular walls seen in the HST image 



observed just 2 years earlier (Nl and SI in Figure 2 of iDoyle et al. 



properties of this jet are in good agreement with those predicted by 



2000h. The qualitati ve 



Livio fc Sokei] (I2OOII ). 
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Motivated by the discovery of the ionized jets in M 2-9, we launched a search for 
coUimated ionized winds in butterfl y nebulae using the VLA a nd the ATCA. We selected a 
total of 16 objects from Table 1 of ISoker fc RappaportI (|2000| ) to search for optically thick 
ra dio cores - a charact e ristic signature of ionized jets. The selected objects were classified 
by ISoker fc RappaportI (120001 ) to have either narrow waists (10 objects) or possibly narrow 
waists (6 objects), whose morphol ogy cannot be deter r nined with absolute certainty based 
on the optical images obtained by ICorradi fc Schwarzl (119951 ) . We are interested in finding 
any central source that is optically thick at short wavelengths where the nebular material is 
already optically thin. Our goal is to identify suitable candidates for follow-up observations 
with higher sensitivity and resolution to reveal the structure of the source. We observe each 
nebula with two wavelengths to calculate its spectral index, which tells us if the emission is 
optically thick. In this paper, we report the results of our search for the radio cores. The 
observations and data reduction are described in Section [2J In Section [31 we report that 
radio emission were detected in 9 sources, and that 5 of those detected have optically thick 
radio cores. Section H] discusses the implication of the high mass loss rates derived from the 
flux densities of the radio cores. Lastly, a summary is given in Section [51 



2. Observations and Data Reduction 

The survey sample comprises 11 northern objects (Declination > —40° for the VLA 
observation limits), IRAS 07131-0147, M 1-16, NGC 2818, NGC 6302, 19W32, HB 5, 
NGC 6537, M 3-28, M 1-91, M 2-48, and NGC 7026, observed with the VLA at 1.3 cm 
and 0.7 cm, and 5 southern objects. He 2-25, He 2-36, He 2-84, Th 2-B, and Mz 3, ob- 
served with the ATCA at 6 cm and 3.6 cm. The observing wavelengths are the shortest 
wavelengths feasible at the time of observations, which give the highest angular resolutions 
possible. The optical properties of the objects are listed in Table [H with columns 1 and 
2 giving the common name and the PN G name based on galactic coordinates, columns 3 
and 4 giving Right Ascension and Declination, column 5 indicating if they appear to have a 
very narrow waist, column 6 stating if the PN has a compact central source visible in optical 
images, and columns 7 indicating if the nebula shows point o r mirror symmetric structure. 



These properties are mainly from ISoker fc Rappaportf (l2000l ). with some update obtained 



from recent optical images. A question mark means the attribute cannot be determined 
from available images. 



- 5 - 



2.1. VLA Observations 

The 11 northern objects were observed at 1.3 cm and 0.7 cm with the A-configuration in 
2002 March 10. The A-configuration is the most extended configuration, hence it provide the 
highest angular resolution achieved with the VLA. It allows us to not just look for but resolve 
if possible a central compact and optically think source, as well as to better discriminate 
against if not resolve out extended nebular emission. The observations were conducted in 
fast-switching mode with a calibration cycle time of just 2-3 minutes in order to correct for 
rapid phase variation at such short wavelengths. The integration time for the individual 
sources is about 20 minutes each at 1.3 cm and 0.7 cm. 3C286 was observed for about 5 
minutes in order to calibrate the fiux scales for individual sources. The data were reduced 
and analyzed using the package AIPS following the standard procedures fo r reducing high- 
frequency data outlined in Appendix D of the AIPS Cookbook (|NRAO||2003| ) . The calibrated 
data were then Fourier transformed, weighted and CLEANed to generate the final CLEAN 
images. Natural weighting was applied in order to obtain images with minimum-noise levels. 
We achieved an angular resolution of ~ 0.08" at 1.3 cm and ~ 0.04" at 0.7 cm for the final 
CLEAN maps. 



2.2. ATCA Observations 

The 5 southern objects were observed at 6 cm and 3.6 cm with the 6-km "6C" con- 
figuration in 2003 March 15 and 16. These two wavelengths were the shortest wavelengths 
feasible when the observations took place. With 6-km being the most extended configura- 
tion, they provide the highest angular resolution achievable at the time of observation. We 
observed each source for 15 minutes roughly every 40 minutes over a period of 12 hours in 
order to optimize the ^i^^coverage, obtaining a total integration time of 4-5 hours on each 
object. The primary calibrator 1934-638 was observed for fiux calibration. The data were 
reduced using the reduction packag e MIRIAD following the procedures outlined in Miriad 



3n[ 



Users Guide ( ISault fc Killeedl2003l ). The final CLEAN images were obtained by perform- 
ing inverse Fourier transforms from the calibrated visibility datasets to images, and then 
deconvolve the sidelobes using the CLEAN algorithm. Natural weighting was used and we 
attained an angular resolution of ~ 3" at 6 cm and ~ 1.5" at 3.6 cm. Where applicable, the 
self-calibration was performed to improve the image quality. In the end, only He 2-84 has 
strong enough emission to allow phase self-calibration. 
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2.3. Optical Observations 

Three objects, He 2-36, He 2-84, and 19W32, which have radio emission detected in our 
study but have no high resolution optical images, were observed with the 2.3m Telescope 
in Siding Spring Observatory (SSO) in 2003 May 26-30. The Imager instrument and three 
narrow-band filters. Ha, [N II] and [O III] were used to take photometric images. The 
center wavelength and bandwidth of each filter are listed in Table [21 along with the best 
achieved FWHM during the observations. These filters' bandwidths of < 10 A allow minimize 
contaminations from adjacent emission lines. The CCD has dimensions of 1024x1024, with 
a focal plane scale of 0.59"/pixel determined by the telescope optics. 



The images were reduced with standard techniques using IRAF software (lMasseylll997l ). 
The averaged bias frames and the overscan level were subtracted from each target image, 
which was then trimmed to discard the overscan region and the first and last few columns 
and rows. The dark current was negligible during the observations, so no dark current 
subtraction was performed. Finally, the images were divided by the flat-field image for each 
fllter. The exposure frames of each object were then combined to generate the final images 
in each filter, using a statistical method that removes cosmic ray contamination. 



3. Results and Analyses 

Among the 16 objects observed, we detected radio emission in nine of them, and five 
of which have a central compact source. In Table [3] we list the detect limits {3a) of our 
observations, and the type of emission that were detected. A more detailed description of 
the detections and the analysis for the compact sources are given in the following subsections. 



3.1. Northern Survey 



Of the 11 objects observed by VIA, we found that NGC 6302, Hubble 5 and NGC 6537 
only exhibit extended radio emission at 1.3 cm that could not be properly mapped due to 
the lack of short baselines in our observations. Two objects, 19W32 and M 1-91, show only a 
compact central source at both 1.3 cm and 0.7 cm. No emission was detected for the remain- 
ing six objects. This is the first time the centers of all these nebulae have been examined 
at such high angular resolutions at radio wavelengths. Note that for the six objects with 
non -detection in our observa tions, some have been detected at 6 cm by IZijlstra et al.l (119891 ) 
and lAaquist fc Kwok Jl990h . whose surveys are more sensitive to the extended emission. 
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Here we present the images for the radio cores of 19W32 and M 1-91, along with their 
optical images (Figures [1] and [2]). Both objects exhibit a compact stellar- like optical core. 
We obtained the optical images of 19W32 from SSO 2.3m telescope observations as de- 
scribed above, those of M 1-91 from the lAC (Instituto de Astrofi'sica de Canarias) catalog 
( iManchado et al.lll996l ). All images are oriented with north up and east to the left. 



3.1.1. 19W32 

The Ha and [N II] images of 19W32 show a prominent stellar-like optical core and two 
lobes in opposite directions, extending ~ 10" to the northeast and southwest. Only the Ha 
image is presented in Figure [H No [O III] emission was detected for this nebula. The radio 
continuum images only show a compact core, which coincides with the optical stellar core. 



3.1.2. M 1-91 

The Ha, [N II] and [O III] images of M 1-91 shows an optical core and two lobes in 
opposite directions, extending at least 20" to the northeast and southwest. Since the images 
of the three filters show similar structures, only the Ha image is presented in Figure [2l M 1- 
91's morphology is very similar to M 2-9, except that the brightenings of the lobes being 
point symmetric instead of mirror symmetric. The radio continuum images only show a 
compact core, which coincides with the optical stellar core. The 1.3 cm emission of M 1-91 
shows some extended structures that might align with the point-symmetric structures. This 
needs further observations to confirm. 



3.1.3. Spectral Indices of the Radio Cores 



The spectral index, defined as the slope of the radio spectrum ^j^-^, indicates the 
emission mechanisms and optical depths of the radio emission. Optically thick free-free 
emission has a rising spectrum, optically thin free- free emission has a flat spectrum (~ —0.1), 
and optically thin nonthermal emission has a steeply falling spectrum. In the case of a 
spherically symmetric nebula, the radio emission has a spectral ii idex of 0.6 {Sy oc 



for an isothermal ste llar wind with a constant expansion velocity (jWright fc Barlow! Il975 



Panagia &: Fellilll975l ). For a special case of coUimated io nized winds with constant opening 
angle, the radio emission also has a spectral index of 0.6 (lReynoldslll986l ). 



In order to determine the properties of the radio cores of 19W32 and M 1-91, we cal- 
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culated the spectral indices of their radio emission. To determine the flux densities of the 
radio cores observed in 19W32 and M 1-91, a two-dimensional Gaussian structure was fit 
to the measurements by applying the AIPS tasks IMFIT to the clean map and OMFIT to 
the visibility data. Because of the low S/N ratio in these observations, neither task gives 
an obviously superior fit. Therefore results from both tasks are given. The entire dataset 
were used for fitting at both 1.3 and 0.7 cm. We also fit a more restricted dataset at 0.7 cm 
chosen to have the same sampling range in angular scales as at 1.3 cm to check for a greater 
inclusion of any extended nebula emission in the larger synthesized beam. We found that 
both procedures gave the same results within measurement uncertainties (< la). Thus, in 
Table HI only the results obtained by fitting the datasets over their entire uv range are listed. 
The rising spectra of both 19W32 and M 1-91 suggest that their radio cores are produced 
by optically thick free-free emission. Fits to the measured visibilities show that these cores 
are resolved along one dimension; i.e., the minor axis of the synthesized beam. The size of 
the emitting source will be discussed further in §4.11 



3.2. Southern Survey 

Of the five objects observed by ATCA, He 2-36 shows only extended emission. He 2-25 
and He 2-84 shows only a compact central source, Mz 3 has both, and no emission was 
detected from Th 2-B. Figures [3] to [H] present the radio and optical images of the four objects 
with detections. For the three objects with radio core detections, their optical images all 



show a compact stellar-like cores. The optical images of He 2-25 were obtained from lLee et al. 



( 120061 ) . those of Mz 3 from the HST archive, and those of He 2-36 and He 2-84 from SSO 



2.3m telescope observations. As before, images are oriented with north up, east to the left. 



3.2.1. He 2-25 

The optical images of He 2-25 display a bright central source and a pair of very faint 
bipolar lobes extending toward north and south (Figure [3]). The radio emission of this 
object is only detected at 3.6 cm, which shows a central unresolved source coinciding with 
the optical core. 
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3.2.2. He 2-36 

The optical images of He 2-36 display a bright central source and an equatorial torus 
with limb brightening arms extending clockwise from the torus edges toward the north and 
south. The brightening resembles an S-shape point symmetrical distribution. We presented 
the Ha image in two different scales in order to show the detailed structures (Figure H]). 
The nebula possesses a more moderate degree of bipolarity (lobes are less collimated) than 
other butterfly nebulae observed in this work. The radio images only show emission from 
the equatorial band. Although the radio images show some interesting features, these are 
not relevant to current investigation. 



3.2.3. He 2-84 

The Ha and [N II] images of He 2-84 show an optical core with an equatorial torus. The 
two lobes extending at right angles to this torus toward the northeast and southwest appear 
to have a point symmetrical geometry. The Ha image is displayed in two different scales 
in order to show the detailed structures (Figure E]). The lobes are very faint and cannot be 
seen in the [O III] image. The radio continuum images only show a central unresolved bright 
source, which coincides with the optical core. 



3.2.4. Mz 3 

The optical images of Mz 3 show a bright core, two approximately spherical bipolar 
lobes, and even more extensive filamentary nebulosities on either side of the central star can 
be seen from its Ha and [N II] images. Only the Ha images is presented in Figure [61 The 
total size of the nebula spans over 50" on the sky. The radio images of Mz 3 sh ow emission 



from the central source and the two lobes, similar to the radio images of M 2-9 (iKwok et al. 



19851 ). The radio central source and the two lobes coincide with the optical core and the 



brightest structures of the optical lobes. 



3.2.5. Spectral Indices of the Radio Cores 

He 2-25 has a compact radio core detected only at 3.6 cm. The fact that no detection at 
6 cm has already indicated that this core is optically thick. Using the detection limit at 6 cm, 
a lower limit of spectral index is estimated to be 1.25. However, a more precise evaluation 
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of the spectral index of He 2-25 has to wait for observations at shorter wavelengths. 

He 2-84 and Mz 3 have radio cores detected both at 3.6 cm and 6 cm, suitable for spectral 
index determination. For He 2-84, we applied the task IMFIT to the clean map and UVFIT 
to the visibility data to determine the flux densities of its radio core assuming the radio core 
to have a two-dimensional Gaussian structure. For Mz 3, the flux density was obtained only 
for the bright compact core in the CLEAN map. This core can be seen in Figure [61 The 
total flux density was obtained by summing up the flux densities of the CLEAN components 
in the radio core region, since the convolution of the Clean components with the CLEAN 
beam might result in the confusion of the core emission and nebular emission. The results 
are listed in Table [5] with la uncertainties. Even with only clean components in the core 
region considered, it seems only part of the diffuse emission is excluded. An analysis of this 
is given later. 

The radio core of He 2-84 has a spectral index significantly more negative than the value 
of —0.1 expected for optically thin free-free emission. This suggests that the emission of the 
observed radio core is dominated by more extended emission that is, as expected, better 
resolved out at shorter than at longer wavelengths. Therefore the extended emission makes 
a stronger contribution at 6 cm than at 3.6 cm, resulting in a spectral index that is even 
more negative than expected for optically thin free-free emission. To illustrate the effects of 
the extended emission on the flux measurements. Figure [7] plots amplitude vs. M?>distance 
for He 2-84 at 6 cm and 3.6 cm. The real component of the amplitudes have been averaged 
in several uv bins. The binning size is wider at longer ui'-distance in order to increase the 
signal to noise ratio at longer baseline. As shown in the figure, the amplitude decreases with 
increasing M?>distance which clear demonstrates that an extended component is present. 

In order to mitigate against contamination from the diffuse nebular emission, we ex- 
cluded the short baselines (< 50 kA) that pick up the extended nebular component. We also 
excluded the longest baselines at 3.6 cm in order to provide the same angular resolution as 
at 6 cm. Because of the noise increases caused by removing uv data, the visibility and map 
fits did not yield reliable results. Consequently, fiux densities were derived by summing the 
fiux densities of the CLEAN components. The results are listed in Table El which indeed 
shows that the extended emission constitutes ~ 90% of the fiux at 6 cm, and ~ 80% of 
the fiux at 3.6 cm for He 2-84. A similar analysis was also performed on Mz 3, since some 
contamination from the extended emission is also possible. The results listed in Tabled show 
that for Mz 3, the fiux contributed from the extended emission is relatively small at 6 cm, 
and there is virtually no contamination at 3.6 cm. 

The differences between the core flux densities and resultant spectral indices obtained 
using the different uv distances suggest that for both He 2-84 and Mz 3, the emission comes 
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from two sources, one compact and one diffuse. In both cases, the more careful analysis 
incorporating the same ^^^^coverage (Table [6]) indicates that an optically thick ra dio core has 



been detected. A separate ATCA investigation of the radio properties of Mz 3 by lBains et al. 



( 120041 ) also concluded that Mz 3 has an optically thick core. 



It is worthwhile mentioning that we did find a compact radio source both in 6 cm and 
3.6 cm inside the primary beam of the Th 2-B observation as shown in Figure [H However, 
after carefully examining the position of this source with the known objects in the field using 
simbadQ, we determined that this is a background radio source not related to Th 2-B, since 
it is about 3' away from the optical postion of Th 2-B. This background radio source has a 
spectral index of —1, characteristic of nonthermal radio emission. It is probably a quasar. 



Discussion 



4.1. The Mass Loss Rate 



The VLA and ATCA radio flux densities obtained in the previous se ctions can be used 
to est ir nate the mass loss rates of the nebulae, following the procedure in IWright fc Barlow 
( 119751 ). IWright fc Barlowl gave the mass loss rate for a spherical wind with a spectral index 
of 0.6 as 



M = 0.095- 



M© yr- 



where v is the velocity of the mass loss flow in km s~^, S*^ is the flux observed at frequency 
u in Jy, u is in Hz, D is the distance in kpc, Z is the nuclear charge, /ig is the mean atomic 
weight per electron, Xg is the ratio of electron density to the number density of ionized H, 
i.e., TT-e = XeUp, and g is the Gaunt factor, which can be approximated in the radio regime as 



TT 



17.7 + In 



3/2 



(z//Hz)Z 



(2) 



where Tg is the electron temperature ( lOsterbrockl Il989l . p. 95) 



The constants Xe and fig can be found by assuming the He to H number ratio y and 
the fraction of He that is singly ionized y', since the contributions from heavy elements are 
negligible due to their low abundance. The singly ionized He to H ratio is y'y and the doubly 
ionized He to H ratio is (1 — y')y if there is no neutral He. The ratio of electron density to 



^http:/ /simbad. u-strasbg.fr/ 
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the number density of ionized H becomes Xe 
weight per electron is 



1 + yy' + 2y{l — y'), and the mean atomic 



l + 4y 



1 + yy' + 2y{l - y') 

For y = 0.11 and y' = 0.5, and fi^ have values of 1.165 and 1.236, respectively. 



(3) 



Assuming Tg = 40 00 K (the temperatu re of free electrons are mainly contributed by 
ionization of hydrogen, iReid fc Menteru 119971 ) and Z = 1, we calculate the mass loss rate as 
a function of distance and wind velocity using the flux densities found in our observations. 
Due to the large uncertainties of the mass loss rate parameters, the differences from the flux 
densities at different wavelengths for the same object are nebligible. Therefore we give a 
best estimate of the mass loss rate for each object in Table [71 

To calculate the actual mass loss rates, the wind velocities and the distances to the 
nebulae are needed. Assuming a typical wind velocity of 1000 km s~^, and using the distances 
available from the literature as given in Tabled we obtain M ~ 7x lO"*' M© yr"^ for He 2-25, 
M ~ 10-5 Mq yr"^ for He 2-84, M ~ 3 x 10"^ M© yr"^ for 19W32, M ~ 9 x 10"^ yr^^ 
for M 1-91 and M ~ 2 x 10^^ Mf;^ yr"^ for Mz 3, r espectively. The mass loss rate for Mz 3 is 
consistent with that derived by lBains et al.l ( 2004J ). They found M ~ 5 — 30 x 10^^ Mq yr^^ 
using the range of D available in the literature. 



We also used the Eq. (11) in lWright fc Barlowl ( 1l975l ) to calculate the size of the emitting 
source, under the spherical wind assumption. The size equation, when incorporated Eq. ([1]), 
reduced to 



r z/ 



cm 



5.83 X 10^' Sy^T-'/^u-'D 
3.9 X lO^^y^T-^/V-i arcsec 



(4) 



where r corresponds to the radius of the emitting region where the optical depth T,y = 0.244. 
For the unresolved sources, Eq. (jl]) predicts a source size smaller than the synthesized beam. 
For the two resolved sources 19W32 and M 1-91, we listed the measured dimensions along 
the resolved axis and the predicted source sizes in Table [SI The predicted sizes are about 
4 or 5 times smaller than the measured dimensions. It is possible that we underestimated 
the predicted size by giving the wrong assumption of temperature in Eq. (jl]). However, to 
make the predicted size consistent with the measured dimension, the temperature has to 
be unrealistically small (~ 200 K), in which case the winds are not ionized. Therefore, we 
concluded that the predicted sizes are smaller is because these sources are elongated, possibly 
by coUimated winds (jets). 



If we approximate a coUimated source as an ellipse with major and minor axes of a and 
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b, and the coUimated wind has an emission area as the same as a spherical wind, then 

a = r^/2/fo. (5) 

where = 26/a is the opening angle of the jet in radians. Using the measured dimension of 
the resolved axis as a, we can derive the opening angle of the collimated jet for 19W32 and 
M 1-91 as listed in Table El The op ening angle can then be used to revise the mass loss rate. 
As pointed out by iReynoldd (119861 ). a collimated wind with a much lower mass loss rate can 
easily produce the same radio intensity as an uncollimated wind with a much higher mass 
loss rate, with the relationship for a spectral index of 0.6 is 



M(jet) 



0.20^n(sinzl 



-1/4 



(6) 



M (spherical wind) 

where i is the inclination of the major axis of the wind to our line of sight. Assuming i = 90°, 
which gives a lower limit of the ratio, we derived the ratio is ~ 0.02 in both cases. Thus we 
revised the mass loss rates of 19W32 and M 1-91 to M ~ 10"^ Mq yr~^ for collimated winds, 
which are about 50 times lower than assuming spherical symmetry. From the analysis of 
central-star wind P-Cygni profiles observed with the International Ultraviolet Explo rer, the 
mass loss rates from central stars of PNs are 10~^ — 10^"^ Mq yr~^ (jPerinottol Il983l ) . These 
values are consi stent with theoretical v alues expected from winds driven by radiation pressure 
on ionized gas (jPauldrach et al.ll2003l ). The mass loss rates we inferred are at least an order 
of magnitude higher than expected. Such high mass loss rates can be explained if the central 
stars of these nebulae are symbiotic star s, as in the case sus pected for M 2-9 whose i n ferred 
mass loss rate is ~ 6 x 10"^ Mq yr"^ jlim fc KwomI2003I ). Similarly, Isains et al.l kooi ) 
also favored a binary companion scenario instead of the superwind of a single progenitor to 
explain the observational properties of Mz 3. 



4.2. PNs or Symbiotic Stars? 



Symbiotic stars are interacting binary systems, typically consisting of a white dwarf 
that accretes material from the wind of a cool red giant companion, sometimes forming a 
disk around the white dwarf. The radiation from the hot star partially ionizes the wind from 
the cool star. A spectrum of the two together thus shows a combination of absorption lines 



from the surf a ce of the c ool re d gia nt and emission lines 



Merrilll (ll919l l. IPlaskettI (ll928l l. and Merrill fc Hummason 



of spectra, these objects were named symbiotic stars by iMerrilll (119581 ). 



r om th e nebula. First noticed by 
(ll932l ) as st ars with combination 



T he poss ibility of some well-k nown PNs can in fact be symbiotic stars was discussed by 



Kwokl (l2003bl ) and ICorradil (120031 ). With the discovery of increasing number of PNs with 
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binary nuclei, the simplistic notion that PNs have single central stars and symbiotics are 
binaries is no longer viable. It is often difficult to distinguish the two classes of objects 
by morphology, spectrum, or spectral energy distribution. Only in cases where periodic 
photometric variability as a result of pulsation, or the observations of photospheric molecular 
absorption band in the spectrum, that one can be certain that we are dealing with a symbiotic 
system. However, it is difficult to rule out that certain PNs, including those under study 
here, are NOT symbiotic stars. Even the evolutionary status of extremely well-studied 
objects such as M 2-9 is still not certain. 

If these nebulae with radio cores are indeed symbiotic star systems, the ionized wind we 
measured would originate from the white dwarf companion. The mass loss rate of this wind 
does not directly tell the mass loss rate of the AGB star, but they must be related since the 
accretion disk is material captured from the AGB star wind. If there is a balance between 
accretion onto the disk and ejection to form a jet, then the mass loss rate of the AGB star 
must be at least as high as the mass loss rate of the jet. Assuming a spherically symmetric 
outflow and only a fraction of the slow neutral wind from the AGB star is captured by the 
WD dwarf, the corresponding mass loss rate of the AGB star would be at least in excess 
of 10~^ Mq yr~^. Such values are seldom observed in AGB stars. This problem can be 
overcome if the AGB star wind is strongly concentrated in the binary orbital plane, as is 
predicted for symbiotic star systems with relatively small orbital separations. 

The remaining objects have no detected radio core emission because they are either too 
distant or too faint. None of those with known distances (0.5 - 7 kpc) are farther than the 
most distant PN with a detected radio core (M 1-91, 7 kpc). If all these objects have the 
same radio luminosity, we should be able to detect all of them. Therefore the more likely 
explanation is that these objects have intrinsically faint radio emission. It is possible that 
those objects evolved from single stars thus do not produce strong radio emission. Since the 
mass loss rate from the central stars of PNs are at least an order of magnitude smaller than 
that of the symbiotic systems, even if these nebulae do have ionized jets from the central 
stars, their radio emission would be too faint to be detected. Thus radio observations may 
provide a way to separate symbiotic nebulae from genuine planetary nebulae. This would 
also mean that there is more than one mechanism to give a nebula a narrow-waist shape. 
For s ingle stars, one of the mor e possible mechanisms for shaping nebulae is magnetic fields 



le.g. 



Garcia-Segura et al.lll999l ). 
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5. Summciry 

We have used the VLA and ATCA to search for signatures of ionized winds - optically 
thick radio cores - in sixteen narrow-waist bipolar nebulae. Eleven nebulae were observed 
with the VLA at 1.3 cm and 0.7 cm with resolutions of 0.08" and 0.04", respectively. This is 
the first time that the centers of these objects have been systematically studied at such high 
angular resolution at radio wavelengths. Two objects, 19W32 and M 1-91, were found to 
exhibit a compact central source with a positive spectral index. Three more objects examined 
with the VLA, NGC 6302, Hubble 5, and NGC 6537, exhibit extended radio emission that 
cannot be mapped properly at such short wavelengths. No emission was detected for six 
other objects in the VLA survey. 

Five nebulae were observed with the ATCA at 6 cm and 3.6 cm with resolutions of 3" 
and 1.5", respectively. He 2-25 was found to have a compact central source only at 3.6 cm. 
He 2-84 and Mz 3 were found to contain a compact central source at both wavelengths. 
When short baselines that are sensitive to extended nebular emission were excluded, we 
found that the spectral index of these two radio cores changes with ■uv-coverage, suggesting 
that an optically thick core has been detected. In addition to its compact radio core, Mz 3 
also exhibits extended nebulae emission that coincides with the two lobes seen in the HST 
optical images. He 2-36 also shows extended nebular emission but no compact core. No 
radio emission was detected for Th 2-B in the ATCA study. 

For 19W32 and M 1-91, we found that the predicted source sizes are smaller than the 
measured dimensions, which we attribute to the radio cores being elongated, implying a 
coUimated wind or jet. The possibihty that these radio cores represent COFs responsible 
for shaping the bipolar morphology of the nebulae is particularly intriguing. The derived 
mass loss rates of the wind from the radio observations are high in comparison to typical fast 
winds from central stars of planetary nebulae, and one possible solution is that the winds 
are not spherical symmetric but instead are collimated. If these winds are collimated by 
a binary system, then the radio imaging offers the best way to test such models. The five 
objects with radio core detections all exhibit a compact stellar-like optical core. This needs 
further investigation in order to understand the connection between the optical morphology 
and the radio ionized jets. 

We are planning follow-up observations of these radio cores with higher angular resolu- 
tion in order to study the structures of the cores. We have recently obtained observations 
of M 1-91 with the VLA A-hPT (Pie Town) configuration. The addition of the Pie Town 
antenna will increase the resolution of the VLA by a factor of two while keeping the full 
sensitivity of the VLA. The new data will add to our knowledge of the core properties of 
M 1-91 once the analysis is completed. 
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Table 1. The properties of the selected PNs with the first 11 objects observed by the VLA 

and the last 5 objects observed by the ATCA. 



Object 


PN G 


RA 


Dec 


VNW 


Optical core 


Point/M: 


irlAo u 1 io i-u 14: 1 




U 1 .ioAZ.o 


— Ui.0z.4z 


yes 


yes 


7 


i\/r 1 1 
M i-io 


ZZO. I + Uo.o 


U 1 .61 .io.o 


— Uy.oo.4o 




7 


yes 


iNvjrU ZOio 


zoi.y + Uo.o 


Uy. io.Ui.o 


—60.6 1 .61 




no 


7 


iNLtU OdUz 


o4i).o + Ui.U 


1 7- 1 Q • /I /I O 

i 1 .16.4:4:. Z 


—6 / .U0.i4 


yes 


no 


no 


1 OAST'QO 

iy Woz 


ooy.z -|- Ui.z 


i / :oy:Uz.o 


— zo:od:oo 


yes 


yes 


7 


HB 5 


359.3 - 00.9 


17:47:56.8 


-29:59:53 


yes 


yes 


yes 


NGC 6537 


010.1 + 00.7 


18:05:13.2 


-19:50:34 


yes 


? 


yes 


M 3-28 


021.8 - 00.4 


18:32:41.3 


-10:06:05 


yes 


? 


yes 


M 1-91 


061.3 + 03.6 


19:32:57.3 


+26:52:43 


yes 


yes 


yes 


M 2-48 


062.4- 00.2 


19:50:27.9 


+25:54:28 


yes 


? 


no 


NGC 7026 


089.0 + 00.3 


21:06:18.5 


+47:51:08 


? 


? 


yes 


He 2-25 


275.2 - 03.7 


09:18:01.3 


-54:39:29 


yes 


yes 


? 


He 2-36 


279.6 - 03.1 


09:43:25.6 


-57:16:56 


? 


yes 


yes 


He 2-84 


300.4 - 00.9 


12:28:46.7 


-63:44:35 


yes 


yes 


yes 


Th 2-B 




13:28:38.2 


-63:49:42 


? 


yes 


? 


Mz 3 


331.7-01.0 


16:17:12.6 


-51:59:08 


yes 


yes 


yes 



Table 2. Optical observing parameters 



Filter Name Center A Bandwidth AA FWHM 

(A) (A) (") 



Ha 6568 8 2.2 

[Nil] 6588 8.8 2.4 

[OIII] 5007 10 2.4 
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Table 3. Results of the VLA and ATCA Surveys. 





Detection limit 


Synthesized beam size 


-L VJJt. *Ji J— /Ctt-L^tl^Jil 




(Jy) 


( 


;") 






1.3 cm 


0.7 cm 


1.3 cm 


0.7 cm 




IRAS 07131-0147 


3.0E-4 


7.2E-4 


0.11 X 0.079 


0.080 X 0.041 




M 1-16 


3.3E-4 


7.8E-4 


0.13 X 0.077 


0.080 X 0.042 




NGC 2818 


5.1E-4 


1.5E-3 


0.26 X 0.082 


0.16 X 0.042 




NGC 6302 


6.6E-4 


4.2E-3 


0.26 X 0.083 


0.17 X 0.041 


extended 


19W32 


3.6E-4 


9.9E-4 


0.19 X 0.082 


0.12 X 0.040 


compact 


HB 5 


3.9E-4 


9.6E-4 


0.20 X 0.078 


0.12 X 0.042 


extended 


NGC 6537 


3.9E-4 


9.9E-4 


0.16 X 0.077 


0.094 X 0.043 


extended 


M 3-28 


2.9E-4 


7.2E-4 


0.14 X 0.079 


0.080 X 0.046 




M 1-91 


3.0E-4 


8.1E-4 


0.10 X 0.083 


0.061 X 0.045 


compact 


M 2-48 


3.0E-4 


7.5E-4 


0.10 X 0.084 


0.060 X 0.045 




NGC 7026 


3.3E-4 


8.1E-4 


0.096 X 0.091 


0.068 X 0.042 






6 cm 


3.6 cm 


6 cm 


3.6 cm 




He 2-25 


3.3E-4 


1.9E-4 


3.74 X 2.71 


1.98 X 1.42 


compact 


He 2-36 


1.8E-4 


1.8E-4 


3.59 X 2.87 


1.88 X 1.52 


extended 


He 2-84 


1.7E-4 


2.1E-4 


3.45 X 3.02 


1.76 X 1.60 


compact 


Th 2-B 


1.7E-4 


2.1E-4 


3.47 X 2.71 


1.82 X 1.42 




Mz 3 


4.5E-4 


4.5E-4 


4.09 X 2.11 


2.24 X 1.16 


compact -|- extended 



Table 4. Flux densities for compact cores of 19W32 and M 1-91 



Object 


Fitting Process 


1.3 cm 
(mJy) 


0.7 cm 
(mJy) 


Spectral Index 


19W32 


Visibility fit 
Map fit 


3.58 ± 0.18 
3.48 ± 0.17 


6.18 ± 0.69 
4.76 ± 0.57 


0.81 ± 0.18 
0.47 ± 0.19 


M 1-91 


Visibility fit 
Map fit 


2.97 ± 0.16 
2.65 ± 0.15 


5.11 ± 0.39 
3.97 ± 0.38 


0.79 ± 0.14 
0.64 ± 0.17 
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Tabic 5. Fhix densities for compact cores of He 2-84 and Mz 3 



Object 


Fitting Process 


6 cm 


3.6 cm 


Spectral Index 






(mJy) 


(mJy) 




He 2-25 


Map pixel sum 


< 0.33 


0.69 ± 0.06 


> +1.25 


He 2-84 


Visibility fit 


12.1 ± 0.16 


9.2 ± 0.16 


-0.47 ± 0.04 


Map fit 


8.6 ± 0.06 


7.7 ± 0.08 


-0.18 ± 0.02 



Mz 3 Map pixel sum 15.3 ± 0.15 18.3 ± 0.15 +0.35 ± 0.02 



Table 6. Compact core fluxes vs. uv-range for He 2-84 and Mz 3 



Object 


uv-range 
(kA) 


6 cm 
(mJy) 


3.6 cm 
(mJy) 


Spectral Index 




50-100 


1.83 ± 0.14 


2.60 ± 0.10 


+0.60 ± 0.15 


He 2-84 


65-100 


1.39 ± 0.16 


1.89 ± 0.12 


+0.52 ± 0.22 




75-100 


1.08 ± 0.18 


1.42 ± 0.15 


+0.47 ± 0.34 


Mz 3 


45-100 
55-100 


13.3 ± 0.11 
10.8 ± 0.14 


19.5 ± 0.11 
20.3 ± 0.11 


+0.65 ± 0.02 
+1.07 ± 0.02 
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Table 7. Mass loss rates and distances of the five objects with radio core detection 



Object 




Distance 


D. Reference 




(Moyr-i) 


(kpc) 




19W32 


~ 5 X 10^9 


3.0 


1 


M 1-91 


~ 5 X 10-9 


7.0 


2 


He 2-25 


~ 2 X 10-9 


2.1 


3 


He 2-84 


~ 5 X 10-9 


1.6 


2 


Mz 3 


~ 3 X 10-^ 


3.3 


4 



tM(M0yr-i) = M^d x (t;/km s"^)(L'/kpc; 



3/2 



References. 



1. 



van der Veen et ahl (119891 ) 



Kohoutekl (Il982l ): 2. iMaciell mm : 3. IPhiUipd (l2004f ): 4 



Table 8. The measured dimensions and predicted sizes of 19W32 and M 1-91 



Object 


A 


Resolved Axis 




riy) 


^0 




(cm) 


(mas) 


n 


(mas) 


n 




1.3 


63 ± 11 


51 ± 19 


17 


7.8 


19W32 


0.7 


51 ± 10 


56 ± 46 


10 


4.3 




1.3 


77 ± 7 


86 ± 6 


14 


4.0 


M 1-91 




0.7 


34 ± 6 


58 ± 12 


9 


8.0 



Positon angle of the resolved axis. 
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17 39 02.38 02.37 02.86 02.35 02.84 02.33 02.82 17 39 02.38 02.37 02.36 02.85 02.84 02.83 02.82 

RIGHT ASCENSION (J2000) RIGHT ASCENSION (J2000) 

Grey scale tlux range= -0.515 2.917 MIIIIJY/BEAM Grey scale tlux range= -1.541 3.209 MIIIMY/BEAM 

Cent peak flux = 2.9166E-03 JY/BEAM Cont peak (lux = 3.2095E-03 JY/BEAM 

Levs = 1.200E-a4 ■ (3, 5, 10, 15, 20, 23) Levs = 3.300E-04 * (3, 5, 9) 



Fig. 1. — Optical and radio images of 19W32. a) SSO Ha image in a linear scale, b) 
VLA 1.3 cm (left) and 0.7 cm (right) radio continuum images with contours superposed. 
The contours start at 3(7, with the noise level cr = 1.1 x 10~^ Jy/beam at 1.3 cm and 
(7 = 3.6 X 10~^ Jy/beam at 0.7 cm. The radio images are on a much smaller scale, with the 
axis range only ~ 1". The position of the radio core coincides with the optical stellar core. 
In this and subsequent figures, the ellipse in the box in the lower-left corner represents the 
half-power extent of the synthesized beam. 



10" 




32 57.53 57.62 57.61 57.60 57.59 57.58 
RIGHT ASCENSION (J2000) 
Grey scale flux ranges -0.431 1.965 MIIIIJY/BEAM 
Com peak flux = 1.9655E-03 JY/BEAM 
Levs = 1.00DE-04 * (3, 5, 10, 15, 18) 




I I \ L 

2 57.63 57.62 57.61 57.60 57.59 57.58 57.57 
RIGHT ASCENSION (J2000) 
Grey scale flux range= -1.166 3.431 MIIIIJY/BEAM 
Cont peak flux = 3.4313E-03 JY/BEAM 
Levs = 2.700E-a4 • (3, 5, 8, 12) 



Fig. 2. — Optical and radio images of M 1-91. a) lAC Ha image in a linear scale, b) 
VLA 1.3 cm (left) and 0.7 cm (right) radio continuum image with contours superposed. 
The contours start at 3(7, with the noise level a — 9.95 x 10~^ Jy/beam at 1.3 cm and 
a — 2.7 X 10~^ Jy/beam at 0.7 cm. The radio images are on a much smaller scale, with the 
axis range only ~ 1". The position of the radio core coincides with the optical stellar core. 
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-54 39 20 
22 
24 




' I \ \ 

09 18 02.0 01.5 01.0 

RIGHT ASCENSION (J2000) 
Grey scale flux ranges -0.254 0.858 MilliJY/BEAM 
Cent peak flux = 8.5751 E-04 JY/BEAM 
Levs = 6.300E-05 ' (3, 5, 7, 10, 12) 



Fig. 3. — Optical and radio images of He 2-25. a) NTT R-band image (iLee et al.ll2006l ) in 
a logarithmic intensity scale, b) ATCA 3.6 cm (right) radio continuum image in greyscale 
with contours superposed. The contours start at 3a with the noise level a = 6.3 x 10~^ 
Jy/beam at 3.6 cm. No emission is detected at 6 cm. 
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09 43 29 28 27 26 25 24 23 22 09 43 29 28 27 26 25 24 23 22 

RIGHT ASCENSION (J2000) RIGHT ASCENSION (J2000) 

Grey scale flux ranges -0.28 10.24 MilliJY/BEAM Grey scale flux ranges -0.287 3.698 MilliJY/BEAM 

Cant peak flux = 1.0244E-02 JY/BEAM Cont peak fJux = 3.6957E-03 JY/BEAM 

Uvs = 6.000E-05 * (3, 5, 10, 20, 50, 80, 120, 160) Levs = 6.100E-05 * (3, 5, 10, 20, 40, 55) 



Fig. 4. — Optical and radio images of He 2-36. a) SSO Ha images show the detailed 
structures of the center in a linear scale (left) and the outer lobes in a logarithmic intensity 
scale (right), b) ATCA 6 cm (left) and 3.6 cm (right) radio continuum image in greyscale 
with contours superposed. The contours start at 3a with the noise level cr = 6.0 x 10"^ 
Jy/beam at 6 cm and cr = 6.1 x 10~^ Jy/beam at 3.6 cm. 
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Fig. 5. — Optical and radio images of He 2-84. a) SSO Ha images show the detailed 

structures of the center in a linear scale (left) and the outer lobes in a logarithmic intensity 
scale (right), b) ATCA 6 cm (left) and 3.6 cm (right) radio continuum image in greyscale 
with contours superposed. The contours start at So" with the noise level a = 5.7 x 10~^ 
Jy/beam at 6 cm and a = 7.0 x 10^^ Jy/beam at 3.6 cm. The position of the radio source 
coincides with the optical stellar core. 
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Fig. 6. — Optical and radio images of Mz 3. a) HST Ha image in logarithmic scales, b) 
ATCA 6 cm (left) and 3.6 cm (right) radio continuum image in greyscale with contours 
superposed. The contours start at 3cr with the noise level cr = 1.5 x 10^"^ Jy/beam at 6 cm 
and (7 = 1.5 X 10~^ Jy/beam at 3.6 cm. The radio core and bipolar lobes coincide with the 
optical stellar core and the brightest structures of the optical lobes. 
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Fig. 7. — Amplitude vs. 'u?>distance for He 2-84 at 6 cm (left) and 3.6 cm (right). The real 
component of the ampplitudes are averaged in several uv bins. The binning size is wider in 
the longer ?/t»-distance in order to increase the signal to noise ratio at longer baselines. 
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Fig. 8. — Images of a background radio source near Th 2-B at 6 cm (left) and 3.6 cm (right). 



